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Abstract 

Nucleon-nucleon potentials evolved to low momentum, which show great promise in 
few- and many-body calculations, have generally been formulated with a sharp cutoff 
on relative momenta. However, a sharp cutoff has technical disadvantages and can 
cause convergence problems at the 10-100 keV level in the deuteron and triton. This 
motivates using smooth momentum-space regulators as an alternative. We generate 
low-momentum interactions with smooth cutoffs both through energy-independent 
renormalization group methods and using a multi-step process based on the Bloch- 
Horowitz approach. We find greatly improved convergence for calculations of the 
deuteron and triton binding energies in a harmonic oscillator basis compared to 
results with a sharp cutoff. Even a slight evolution of chiral effective field theory 
interactions to lower momenta is beneficial. The renormalization group preserves 
the long-range part of the interaction, and consequently the renormalization of long- 
range operators, such as the quadrupole moment, the radius and (1/t), is small. This 
demonstrates that low-energy observables in the deuteron are reproduced without 
short-range correlations in the wave function. 



1 Introduction 

Internucleon potentials with variable momentum cutoffs, known generically as 
"Mowfc," show great promise for few- and many-body calculations [1,2,3,4,5,6,7]. 
Changing the cutoff leaves observables unchanged by construction, but shifts 
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contributions between the potential and the sums over intermediate states 
in loop integrals. These shifts can weaken or largely eliminate sources of 
non-perturbative behavior such as strong short-range repulsion or the ten- 
sor force [8]. An additional bonus is that the corresponding three-nucleon 
interactions become perturbative at lower cutoffs [4]. As a result, it is found 
in practice that few- and many-body calculations can be greatly simplified or 
converge more rapidly by lowering the cutoff. This has been observed with 
few-body variational methods [9,10], the coupled-cluster approach [11], and 
for nuclear matter [5]. 

Low-momentum nucleon-nucleon interactions were originally derived and con- 
structed in energy-independent form using model-space methods (such as Lee- 
Suzuki [1,2] or Okubo [12]) applied in momentum space. These approaches 
define orthogonal subspaces with projection operators P and Q, such that 
P + Q = 1 and PQ = QP = 0. In momentum space, the latter condition 
implies a sharp cutoff A in (relative) momentum, so that P-space integrals 
run from to A, while Q-space integrals run from A to oo (or to a large 
"bare" cutoff). Subsequently, this V\ ow k construction was shown to be equiva- 
lent to a Renormalization Group (RG) treatment, derived by requiring cutoff 
independence of the half- or fully-on-shell T matrix [3]. With a sharp cutoff, 
the equations take a particularly simple form and two-body observables are 
preserved for all momenta up to the cutoff. 

However, a sharp cutoff also leads to cusp-like behavior for the interaction 
(close to the cutoff) in some channels and for the deuteron wave function, 
which becomes increasingly evident as the cutoff is lowered below 2fm _1 . In 
some applications, this leads to slow convergence, for example at the 10- 
100 keV level in few-body calculations using harmonic oscillator bases (see 
Fig. 1). The reduction of the repulsive short-range interaction simultaneously 
reduces short-range correlations in the wave functions, which means that vari- 
ational calculations can be effective with much simpler trial ansatze [9]. One 
would expect that such calculations for the deuteron and the triton, which are 
particularly low-energy bound states, should show improvement for cutoffs well 
below 2fm _1 , but instead a degradation was observed in Ref. [9]. This result 
was attributed to the use of sharp cutoffs and a preliminary study [10] showed 
that these problems are alleviated by using a smooth-cutoff low-momentum 
interaction. In this paper, we verify this conclusion and explore in detail the 
construction and application of V\ ow k interactions with smooth cutoffs. 

While smooth cutoffs seem incompatible with methods requiring PQ = 0, it is 
not a conceptual problem for the RG approach. Indeed, there is an appreciable 
literature on smooth-cutoff regulators for applications of the functional or ex- 
act RG [14]. The functional RG keeps invariant the full generating functional, 
which translates into preserving all matrix elements of the inter-nucleon T 
matrix. While this straightforwardly leads to RG equations, it also implies 
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Fig. 1. The triton binding energy Ef calculated from a direct diagonalization in a 
harmonic oscillator basis of the low-momentum Hamiltonian derived from the Ar- 
gonne v±s potential [13] with cutoff A = 2fm _1 , as a function of the size of the 
oscillator space (A max hu> excitations). The open circles are calculated with a sharp 
cutoff for a fixed oscillator parameter b while the filled ones correspond to optimiz- 
ing b at each A max . The dashed line indicates the exact Faddeev result using the 
sharp-cutoff interaction [4], and shows the slow convergence of the diagonalization 
at the 100 kev level. The squares are for a smooth Fermi-Dirac regulator, Eq. (32), 
that solves the convergence problem. 

an energy-dependent interaction, which is undesirable for practical few- and 
many-body calculations. We resolve this conflict in Sect. 2 by constructing a 
low-momentum, energy-independent interaction with smooth cutoffs in three 
steps: 



(1) Evolve a large-cutoff potential to a lower, smooth cutoff while preserv- 
ing the full off-shell T matrix. This generates an energy- dependent low- 
momentum interaction. 

(2) Convert the energy dependence to momentum dependence, which results 
in a non-hermit ian smooth-cutoff interaction. 

(3) Perform a similarity transformation to hermitize the low-momentum in- 
teraction. 



Along with the possibility of many different functional forms for the smooth 
regulators, various hermitization schemes are possible, which reflects the gen- 
eral freedom in low-energy effective theories [15]. 

The freedom in defining low-energy potentials has consequences in practical 
calculations. This is already evident in Fig. 1, where the (converged) triton 
binding energy for a sharp cutoff is 50keV less than for a particular Fermi- 
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Dirac regulator. The difference reffects the different contribution from the 
(neglected) short-range three-body force. Just as changing a cutoff with a 
sharp regulator moves one along a Tjon line [4], we expect that changing the 
form of the regulator and the hermitization scheme at fixed cutoff also does. 
Thus, each combination of regulator (specified by one or more "sharpness" 
parameters) and hermitization scheme will have different corresponding three- 
body (and higher many-body) interactions. 

In Sect. 3, we derive alternative energy- independent RG equations that are 
generalizations from sharp cutoffs. This approach has the advantage of a one- 
step construction, but accurate solutions of the resulting RG equations for 
low-momentum cutoffs are technically more challenging. The low-momentum 
interactions with smooth cutoffs, with different regulator types and different 
hermitization schemes, are applied in the two-body sector and to calcula- 
tions of the triton in Sect. 4. We test the convergence properties and doc- 
ument the distortions of various combinations. In most cases we derive the 
low-momentum interactions starting from chiral effective field theory (EFT) 
potentials at N 3 LO [16,17], but also use the Argonne v\$ potential [13] for 
comparison. We summarize our conclusions in Sect. 5 and give an outlook for 
future applications. 



2 Smooth cutoff interactions via an energy-dependent RG 

Our goal is to construct a smooth cutoff version of the energy-independent 
and hermitian low-momentum interaction V\ ow k- In this section, we describe a 
three-step method that utilizes an energy-dependent RG equation to lower the 
cutoff, followed by two transformations that remove the energy dependence 
and the resulting non-hermiticity in the interaction. In the next section, we 
derive an equivalent method that uses a hermitian and energy-independent 
RG equation to construct the low-momentum interaction in one step. 

First, we derive the RG equation for an energy- dependent low-momentum 
interaction V e s(E), which is cut off by smooth regulators. It is convenient and 
efficient for numerical calculations to define the partial-wave interaction V c g 
and the corresponding T c s matrix in terms of a reduced potential v and a 
reduced t matrix as 



V cS (k',k;E)=f(k')v(k',k;E)f(k), 
T cS (k',k;E)=f(k')t(k',k;E)f(k), 



(1) 
(2) 
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where f(k) is a smooth cutoff function satisfying 

/(jfe)*^l and /(Jfe)^O. (3) 

The regulator functions are the same in each partial wave and possible choices 
are discussed in Sect. 4. 

Given the low-momentum interaction V e g (E), the reduced fully-off-shell t ma- 
trix satisfies the Lippmann-Schwinger equation, 

t{k,k;E)=v{k,k;E) + - p dp , (4) 

n Jo E — p z 

where we use units with h — c — m — 1 and a principal value integral 
is implicit here and in the following. Note that the smooth cutoff is on the 
loop momentum but not on external momenta, and that v and / are cutoff 
dependent. We impose that the fully-off-shell t(k', k; E) is independent of the 
cutoff, dt(k', k; E)/dA = 0. This leads to an energy- dependent RG equation for 
the change in the reduced low-momentum interaction v(k', k; E) as the cutoff is 
lowered. In operator form, we have t — v (l + G^t) and thus v — t (1 + Gq 
With dU- 1 = -U- 1 dU U~\ we obtain (see also Ref. [18]) 

*=- (( l +G f ( )-.^ t (l + Gf( )- = -.^.. w 



Restoring the momentum and energy arguments, the energy- dependent RG 
equation reads 



d .J. hlA 2 r 2 v(k>,p;E)±[p(p)}v(p,k;E) 

— v{k,k;E) = - p dp ^ . 6 

dA n Jo p 1 — E 



For the implementation, we discretize the momentum space on a set of Gaus- 
sian mesh-points, which results in coupled RG equations that can be nu- 
merically integrated using standard methods. It is convenient to define dis- 
cretized plane-wave states as \k) = k^2wk/n \k), where are the Gauss- 
Legendre weights. The new basis states are normalized as (k\p) = 5^^. Tak- 
ing matrix elements of v(E) between the discretized plane- wave states gives 
Vk',k-,E = f k'ky/Wk'Wkv(k' , k] E), and Eq. (6) takes the simple matrix form 

V k ',k;E ~ 2^ V k',p;E p2 _ £ V PME • {' ) 

Here and in the following, subscripts denote discrete momentum labels. In 
addition, it is convenient to convert the principal value integration to a normal 
integration using the standard principal- value subtraction method [19]. 
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If we take a nucleon-nucleon (NN) potential model Vnn as the large-cutoff 
initial condition and numerically integrate the RG equation, the resulting 
energy-dependent v(k', k; E) preserves the fully-off-shell T NN for all external 
momenta and energies, t(k',k;E) = T NN (k',k; E). Therefore, V e g(k',k;E) = 
f(k') v(k', k; E) f(k) preserves the low momentum fully-off-shell T NN matrix up 
to factors of the smooth cutoff function, T e s(k', k; E) = f(k') T NN (/c', k; E) f(k). 

In the second step, we convert the energy dependence of V e a(E) to momen- 
tum dependence by using a method similar to field redefinitions. For this pur- 
pose, we introduce an energy-independent (but non-hermitian) V\ ovr k(k',k) 
that reproduces the half-on-shell T eS matrix (and hence wave functions) as 

V eS (k',k;E), 

(k'\T cS (p 2 )\p) = (k'\V cS (p 2 )\ Xp } = (k'\V lowk \ Xp ) , (8) 

where \x P ) are the eigenstates of the energy-dependent Hamiltonian H cS (p 2 ) = 
T + V e s(p 2 ) with relative kinetic energy T. Using the completeness of the 
interacting eigenstates, this leads to 

V lowk (k',k)= p 2 dp- p' 2 dp'V cS (k',p';p 2 ) Xp (p')x;(k). (9) 

JO 7T JO 

The energy dependence of H c s(p 2 ) necessitates the use of bi-orthogonal com- 
plement vectors (x P \ in the completeness relation, / p 2 dp \Xp){Xp\ — 1- If bound 
states are present, the integral over the continuous scattering states includes 
a summation over bound states as well. 

For the numerical solution of Eq. (9), we first obtain the eigenstates of the 
energy-dependent H ef[ by self-consistently diagonalizing the discretized eigen- 
value equation 

£(4' fc k 2 + Vf :k . Ev ) Xk ,p = E p X w,p ■ (10) 

k 

Next, the discretized bi-orthogonal complement vectors are obtained from the 
matrix equation J2 P Xk, P Xk', P — 8k,k'- Finally, Eq. (9) is evaluated by simple 
matrix multiplication 

V k?J k = V k? P ';E p Xp',p Xk, P ■ (11) 
P,P' 

This procedure is straightforward in practice. The only subtlety arises when 
the discretization of the defining equations for the bi-orthogonal complement 
vectors is such that some of the \x p ) vectors appear to be linearly depen- 
dent. In this case, singular- value-decomposition methods are helpful to solve 
the singular system of equations by setting to zero the problematic small sin- 
gular values when taking the inverse of the Xk, P matrix [19]. The resulting 
energy-independent Vi owk is non-hermitian, and is identical to the solution 
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of the energy-independent RG equation derived from half-on-shell t matrix 
equivalence in the next section. 

The method can be improved further by integrating the energy-dependent RG 
equation, Eq. (6), formally instead of a numerical integration. That is, using 
dv _1 = — v _1 dvv^ 1 = dG$, we obtain v _1 — V^J = Gq — Gq, and we recover 
the Bio ch- Horowitz equation generalized to a smooth cutoff, 

IT Jo hi — p z 

Converting the RG equation to this integral equation and solving by matrix 
methods speeds up the calculation considerably and reduces the accumulation 
of numerical errors. 1 

Finally, we note that the initial energy-dependent RG equation is not needed if 
one starts directly from the smooth-cutoff generalization of the Bloch- Horowitz 
equation. For this purpose, we separate the free two-nucleon propagator into 
a smooth-cutoff low-momentum Gq and a high- momentum part Gq. The 
Lippmann-Schwinger equation then reads Tnn(-E) = Vnn + Vnn (Gq(E) + 
Gq(E))Tnn(E), and a rearrangement gives the fully-off-shell T NN (E) matrix 
with only low-momentum propagators, 

T NN (£) = v(E) + v(E) G£(E) T nn (£) , (13) 

where v(E) is the solution to the smooth-cutoff Bloch- Horowitz equation, 
Eq. (12). 

2.1 Hermitian low-momentum interactions 

In the third step, we remove the non-hermiticity of Vi ow fc by a similarity 
transformation that orthogonalizes the set of eigenvectors {|Xp)} °f the non_ 
hermitian Hamiltonian H\ owk = T + Vi ow fc. Following Holt et al. [20], we define 
a transformation Z by 

Z \X P ) = \Q and (£p|£ P ') = <W> ( 14 ) 

where the Kronecker delta normalization implies the discretization procedure 
of the previous section has been carried out. The hermitian low- momentum 
interaction is then given by 

V\ wk = Z H\ ow k Z 1 — T. (15) 

1 One recovers the folded diagram series, if one applies the same trick to the energy- 
independent RG equation discussed in the next section. However, this series does 
not sum to a simple linear integral equation. 
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There is not a unique choice for the transformation Z, which is a reflection of 
the general freedom in low-energy effective theories. In Ref . [20] , it was shown 
how several common hermitization methods correspond to different choices for 
Z . For example, within the Lee-Suzuki framework [21,22] for deriving effective 
interactions (which implies a sharp cutoff corresponding to orthogonal PQ = 
projection operators), the eigenstates of the non-hermitian V\ ow k obey 

(X P \P + u j uj\xp>} = 5 pp/ , (16) 

where uj = QuP is the wave operator that parameterizes the Lee-Suzuki 
decoupling transformation. Identifying Z^Z = P + cu^cu defines a class of valid 
transformations for Z . For example, setting Z = y/P + uo^uo corresponds to 
the hermitization procedure of Okubo [23] and Okamoto and Suzuki [24]. 
Alternatively, one can perform a Cholesky factorization of the symmetric and 
positive-definite operator LlJ = P + cu^cu, where L is the lower-triangular 
Cholesky matrix. Then, Z — L' corresponds to the hermitization method of 
Andreozzi [25]. Another method discussed by Holt et al. is a Gram-Schmidt 
orthogonalization to construct the set of directly from {|x P )}. This 

allows Z to be calculated directly from Eq. (14). Even within the Gram- 
Schmidt method, Z is not unique, due to the freedom in choosing the starting 
vector. All hermitization methods result in low-momentum interactions l^i ow fc 
that preserve the low-momentum fully-on-shell Tnn matrix, up to factors of 
the regulator function T lowfc (/c, k; k 2 ) = f 2 (k) T NN (k, k; k 2 ), and the deuteron 
binding energy. The corresponding three-body interactions will differ, however, 
as discussed below. 

The Gram-Schmidt method can be applied directly to the smooth-cutoff V\ ow k- 
As in Ref. [20], the orthogonal basis {|£ p )} is constructed via 

I6> = ^iilxi> 

|6)=^2l|Xl>+^22|X2> (17) 

Ifs) = Zu\Xi) + Z32IX2) + Z 33 \x 3 ) 

where the Z pp i are determined sequentially so that (£ p |£ p ) = S pp >. We have 
chosen to take the eigenstate of i?i ow fe with lowest energy as the starting vector, 
although any other linear combination could have been used. In practice, the 
modified Gram-Schmidt algorithm with re-orthogonalization of Ref. [26] is 
utilized to guard against round-off errors. The transformation Z corresponding 
to the Gram-Schmidt hermitization is then given by 

Z = J2\Q(X P \ and Z- 1 = £ \ Xp ) . (18) 
p p 
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In contrast to the Gram-Schmidt method, the other hermitization schemes are 
formulated in terms of the Lee-Suzuki wave-operator u, which apparently relies 
on the use of orthogonal projection operators PQ = corresponding to sharp 
cutoffs. In order to generalize the Okubo and Andreozzi hermitization schemes 
to smooth cutoffs, it is necessary to eliminate all references to 00. This is easily 
done by noting that the bi-orthogonal complement vectors are defined through 
(Xp\Xp') = $pp'- For a sharp cutoff, Eq. (16) implies \x P ) = (P + u^u)\x P ), and 
thus 

P + Jlu = J2\Xp)(Xp\- (19) 

p 

For smooth cutoffs, the obvious generalization is to construct the operator 
J2 P \Xp) (Xp\ an d decompose it as 

ztz = £ \x P ) (x P \ • (20) 

p 

As for sharp cutoffs, the generalized Okubo transformation corresponds to 
Z = ^Jj2 P \ Xp) (Xpl; where the square root is taken in the eigenbasis of the 
positive-definite operator J2 P \x P ) (x P \- Similarly, the Andreozzi hermitization 
can be obtained from performing the appropriate Cholesky decomposition. 



3 Energy-independent RG using smooth cutoffs 

In this section, we generalize the energy-independent RG equation of Ref. [3] to 
smooth cutoffs. For energy-independent interactions, the reduced half-on-shell 
t matrix obeys a Lippmann-Schwinger equation with loop integrals smoothly 
cut off by f 2 (p), 

t(k ,k;k ) =v(k ,k) + - p dp , 21 

71 Jo k 2 — p 2 

where the energy-independent low-momentum interaction Vj wfc an d the cor- 
responding half-on-shell Ti ow fc matrix are given by 

V lowk (k',k) = f(k')v(k',k)f(k), (22) 
T lowk (k', k; k 2 ) = f(k') t(k', k; k 2 ) f(k) . (23) 

We again note that the smooth cutoff regulator is on the loop momentum but 
not on external momenta. 

Analogous to the derivation for a sharp cutoff [3] , we impose that the reduced 
half-on-shell t matrix is independent of the cutoff, dt(k' ', k; k 2 ) / dA = 0. This 
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choice preserves the on-shell t matrix while also maintaining energy indepen- 
dence. The result is 



2 r°° 

— / p 2 dp 

TT Jo 



dv(k', p) 
dA 



7r5(p-k) f 2 (p)t(p,k;k 2 ) 
2k 2 k 2 -p 2 

2 r v(k>,p)±[f 2 (p)}t(p,k;k 2 ) 



/ p 2 dp ■ 

TT JO 



p2 _ fc2 



(24) 



Next, we express the term in the square brackets in Eq. (24) by the exact scat- 
tering state \xk) of Viowfc- From ViowfelXfe) = Vi owk \k) + V\ owk G (k 2 ) Ti owk \k), it 
follows that (p\Xk) = (p\k) + [f(p) t(p, k; k 2 ) f{k)]/{k 2 — p 2 )- We can therefore 
write Eq. (24) as 



2 r°° 
- / p 2 dp 

TT JO 



dv(k',p) f{p) 
dA f(k) 



2 r°° 
Xk(p) = ~ p 2 dp2 

TT JO 



df(p) v(k',p) f(p)t(p,k;k 2 ) 



dA 



p 2 — k 2 



or equivalently 



2 foo dv{k',p) 2 foo 

- P dp — "tt f(p) Xk{p) = - P dp 2 

7T JO OA TT JO 



df{p) v{k',p) 
dA p 2 -k 2 



\V\owk\Xk) 



(25) 



I r p 2 dp2 d Mv{k' ) p) ( P \V lowk G(p 2 )\xk), 

TT JO OLA 

(26) 



where G(E) — (E — Hi owk ) 1 is the interacting two-nucleon Green's func- 
tion. Using the completeness relation, / k 2 dkx k (p') Xkip) = ^S(p — p') /p 2 and 
V Umk G(E) = T ]owk (E)Go(E),wehave 



-?L v (k',k) = - [ p 2 dp 
dA TT Jo 



2 r^^v( k',p)±[f 2 (p)}t(p,k;p 2 ) 
p 2 — k 2 



(27) 



which describes the evolution of the reduced low-momentum interaction with 
the cutoff. 



If we take an energy-independent NN potential as large-cutoff initial condition 
and numerically integrate the RG equation, then the resulting v preserves the 
half-on-shell Tnn matrix for all external momenta, t(k', k; k 2 ) = Tnn(&', k; k 2 ). 
Therefore, V\ owk (k f , k) = f(k') v(k', k) f(k) preserves the low-momentum half- 
on-shell Tnn matrix up to factors of the smooth cutoff function. In the limit 
f(p) ~^ 9(A—p) and thus d[f(p) 2 ]/dA — > 5(A—p), we recover the RG equation 
for a sharp cutoff [3] 

d„ 2 \/ lowfc (^,A)T lowfc (A,A:;A 2 ) 
dA Vl ™ k{k ' k) = It T^WW • (28) 
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Finally, we can use the Okubo transformation to hermitize V\ owk . In order to 
generalize the hermitization to smooth cutoffs, we consider the sharp-cutoff 
Okubo transformation under an infinitesimal change of the cutoff. In this 
case, one can show that the RG equation for the hermitian V\ owk is given by 
a symmetrized version of Eq. (28), 



aA 7i 



~V lowk {k',A)T lowk {A,k;A 2 ) 



1 - (fc/A) 2 



+ 



J- low k 

(k',A;A 2 )V lowk (A,k) 
1 - (k'/A) 2 



(29) 



A simple generalization of the Okubo transformation to smooth cutoffs is 
therefore obtained by symmetrizing the smooth-cutoff RG equation, Eq. (27), 
to obtain 



aA 7r Jo 



p 2 dp 



p2 _ fc2 
+ 



t(k',p;p 2 )^[f 2 (p)]v(p,ky 

p2 _ 



(30) 



As in the previous section, the hermitian low-momentum interaction \^i ow fc 
preserves the low-momentum fully-on-shell Tnn matrix, up to factors of the 
regulator function T\ owk (k, k; k 2 ) = f 2 (k) T NN (/c, k; k 2 ), and the deuteron bind- 
ing energy. The freedom in the hermitization method can be expressed as an 
auxiliary condition dt{k! , k; k 2 ) / dA = (k 2 — k' 2 ) $(/c' ,/c), where §(k',k) is a 
function with lim k ^ k /(k 2 — k' 2 ) $(/c', k) = 0. The above RG equation derived 
from the Okubo transformation makes a particular choice for $(&/, k). 

The numerical solution of the energy-independent RG equation is compli- 
cated by the t matrix calculation involved in each step. The computational 
overhead slows down the ODE solver significantly. In addition, the RG equa- 
tion involves two-dimensional interpolations and principal- value integrals over 
narrowly peaked functions. Therefore, it is easy to introduce small errors at 
each step that can accumulate as the cutoff is lowered. Finally, some potential 
models exhibit spurious resonances (at order GeV energies and momenta) [27] 
and these need to be subtracted before solving the energy-independent RG 
equation for these potentials. Because of these difficulties, we have exclusively 
used the three-step method from Sect. 2 to generate the results presented in 
the next section. 
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4 Results 



In this section, we apply the formalism discussed in Sect. 2 to derive hermitian, 
low-momentum interactions with smooth cutoffs. Electromagnetic contribu- 
tions are included in the evolution to low momenta. Since all of the following 
results are for the hermitian V\ ow k, we drop the overbar hereafter. In select- 
ing a smooth regulator function satisfying the conditions of Eq. (3), there is 
obviously much freedom, which parallels the freedom of field redefinitions in 
low-energy effective theories and the functional RG [28]. However, there are 
trade-offs in the choice. A smoother cutoff will dampen more the artifacts of 
a theta-function regulator but will distort more the phase shifts for momenta 
near the cutoff. 

We present results for two choices for f(k), each with a range of parameters. 
These are the exponential form used in current chiral EFT potentials [16,17] 
with integer n determining the smoothness, 

/(*) = e~^ n , (31) 

and a Fermi-Dirac form with a sharper cutoff achieved with smaller e, 

1 



/(*) 



1 _|_ e (fc2-A2)/ £ 2 



(32) 



It may be interesting in the future to explore other choices for the regu- 
lator function, but the general features and advantages are covered by the 
above choices. Other possible regulators (used in different applications, e.g., 
see Ref. [14]) include a power law form with integer n, 



/(*) 



1 



1 + (k 2 /A 2 ) n ' 



(33) 



a hyperbolic tangent form with an e parameter that plays a similar role as in 
the Fermi-Dirac function, 



/(*) = 5 



1 + tanh 



'A 2 - k 2 
Ake 



(34) 



a complementary error function with e parameter, 

/(*) = 




(35) 



and a Strutinsky averaging with e parameter, 

/(*) = \ 



erf — + erf — 



(36) 
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Fig. 2. Plots of the exponential and Fermi-Dirac regulators squared as a function of 
momentum k for A = 2fm~ 1 and a range of parameters n and e. 



In each case, the function and its derivatives are continuous, and the parameter 
n or e controls the smoothness. 

In Fig. 2, the quantity [f(k)] 2 is plotted against k for fixed cutoff A = 2fm _1 
using some candidate parameterizations of the exponential and Fermi-Dirac 
forms. While a sharp cutoff (for which f(k) = 8(A — k)) preserves the on-shell 
T matrix up to the cutoff, the T matrix for a smooth cutoff is multiplied by 
lf(k)} 2 , leading to distortions in the phase shifts near the cutoff. The expo- 
nential regulator ("exp") from Eq. (31) with n = 3 corresponds to what is 
used in N 3 LO chiral potentials [16,17]; it is evident that this regulator applied 
in the present context will significantly distort the phase shifts for momenta 
well below A. As n is increased, the regulator gets sharper; for numerical rea- 
sons, n = 10 is probably the practical upper limit. The Fermi-Dirac form 
( "FD" ) interpolates smoothly between sharp and smooth as a function of the 
e parameter. It causes less distortion for lower momenta than the exponential 
regulator. 

The effects of different regulators on V\ ow k potentials with the same starting 
("bare") potential and the same cutoff are illustrated in Fig. 3. In channels 
where the potential is close to zero at A, such as the 3 Si partial wave, the 
differences between sharp and smooth are slight, particularly as the regulator 
gets sharper (n exp is smoother than £fd = 0.5 fm -1 ; see Fig. 2). However, the 
difference in other channels can be striking, as seen in Fig. 3 for the 3 D 2 partial 
wave. The observed cusp-like behavior is due to reproducing the phase shifts 
for momenta up to the cutoff. The existence of sharp regulator artifacts is not 
a problem in principle, as the potential is not an observable, but in practice 
it can lead to convergence problems at the 10-100 keV level in the deuteron 
and triton. 
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One of the striking properties of Vi ow fc with a sharp cutoff is the "collapse" 
of different high-precision NN potentials to almost the same low-momentum 
potential as the cutoff is lowered below 3fm _1 [2]. This behavior is expected 
as a consequence of the same long-range pion-exchange interaction together 
with phase shift equivalence of the potentials up to k ~ 2.1 fm -1 . Therefore, 
it is not surprising that it remains a property of V\ ow k interactions derived 
with smooth regulators, as shown in Fig. 4 for a set of chiral potentials with 
different initial ("bare") cutoffs. Note that the low- momentum interactions 
from different starting potentials are close but not identical. The differences 
will be paralleled by differences in the corresponding short-range three-body 
interactions. 



Low-momentum interactions with smooth cutoffs reproduce the initial phase 
shifts up to factors of the regulator function. The error in phase shifts due 
to the regulator alone is illustrated in Fig. 5 for some representative two- 
body phase shifts as a function of laboratory energy. In particular, for each 
energy, the on-shell T-matrix from a bare potential (in this case the N 3 LO 
chiral potential from Ref. [16]) is multiplied by [f{k)} 2 using the momentum 
k corresponding to that energy. This corresponds to the distortion that would 
be present if numerical errors in constructing V\ ow k were negligible. The latter 
are documented next and are small. We have no universal rule for deciding 
whether a distortion is acceptable; it depends on how it propagates to the 
observable in question. For example, for the low-energy bound-state of the 
deuteron, none of the distortions in Fig. 5 is important. The distortion is 
analogous to the error band from a chiral EFT truncation (but we have not 
formulated a corresponding power counting rule). Therefore, we expect that 
there is no concern if the distortion is comparable to the EFT truncation 
error. In addition, for low-energy properties, the error incurred here can be 
absorbed by the short-range part of the corresponding three-body interactions. 
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Fig. 4. The "collapse" of Vi ow fc interactions derived from N 3 LO chiral potentials 
with a Fermi-Dirac regulator (eFD = 0.5 fm -1 ) as the cutoff is lowered in the 1 So 
and 3 Si channels. The diagonal matrix elements V\ ow k(k,k) are shown, but as in 
Ref. [2], we find similar results for the off-diagonal matrix elements. The different 
lines correspond to different starting potentials with the corresponding cutoffs, A [16] 
or A/A [17], in MeV given in the legends. For Ref. [17], A is the spectral function 
cutoff. 



Consequently, the propagation of phase-shift errors to many-body observables 
can only be studied after including three- nucleon forces. It is evident, however, 
that n exp = 8 and epD = 0.5 only distort minimally. We will show below 
that these are also good choices for convergence. In future applications, the 
regulator effects can be tested by varying the parameters determining the 
smoothness. 
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Fig. 5. The effects of different regulators on the initial phase shifts 5 and the mixing 
parameter e\ as a function of laboratory energy is shown for selected channels. 
The phase shifts are calculated using the N 3 LO chiral potential from Ref. [16]. 
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Fig. 6. Errors in phase shifts relative to regulated bare phase shifts for three hermi- 
tization methods using the exponential regulator with A = 2fm~ 1 . The potentials 
used on the left were derived from Argonne v±s [13] with n ex p = 4 and those used on 
the right were derived from the N 3 LO chiral potential from Ref. [16] with n, 
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Phase shifts can also differ from those calculated from the input ("bare") 
potential because of numerical errors in generating the low-momentum in- 
teraction. In Fig. 6, the effects of numerical errors on the phase shifts for 
different hermitization methods (but fixed regulator) are isolated by plotting 
the difference of the calculated low-momentum phase shifts and the distorted 
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Fig. 7. S-state and D-state components of the deuteron wave function in momentum 
space (u(k) and w(k) respectively) for the bare N 3 LO chiral potential from Ref. [16] 
and those derived using smooth and sharp cutoffs at A = 2.0 fm -1 and 1.5 fm - . 

bare phase shifts. We conclude that the Okubo hermit izat ion is numerically 
more robust than the Gram-Schmidt or Cholesky methods, with errors in the 
phase shifts of order 10~ 4 degrees below 200 MeV and then varying up to 10~ 3 
degrees depending on the regulator and the initial potential. Consequently, 
we use the Okubo hermitization in the following unless otherwise specified. 
For fixed hermitization but different regulators, we have found that relatively 
smooth cutoffs (e.g., epD = 0.5 fm -1 ) achieve 10~ 4 degree accuracy with a 
moderate (of order 50) Gauss points, but sharper cutoffs (e.g., £fd = 0.2 fm^ 1 ) 
have errors for energies above = 200 MeV that can grow as large as 10 _1 
degrees. Greater accuracy can be obtained with a more carefully prescribed 
distribution of points. 2 

The deuteron wave functions for smooth and sharp cutoffs are contrasted in 
momentum space in Fig. 7 and in coordinate space in Fig. 8. We follow the 
notation of Ref. [29], with S-wave and D-wave components denoted in coordi- 
nate space by u and w respectively, and with tildes in momentum space. They 
are normalized as J °° dr [u(r) 2 + w(r) 2 } = 1 and J °° dkk 2 \u(k) 2 + w(k) 2 ] = 1. 
The sharp-cutoff wave functions develop cusp-like structure in momentum 
space below 2fm _1 (see inset in Fig. 7), which are removed by the Fermi- 
Dirac (or any other smooth) regulator. The different momentum-space behav- 
ior is evident in the coordinate-space wave functions as smaller amplitudes 



2 We note that the present scheme occasionally and unsystematically leads to nu- 
merical "glitches" in the low-momentum interaction, particularly for cutoffs close 
to the bare cutoff. These are manifested as discontinuities in the potential and 
are signaled by large discrepancies in calculated matrix elements. We have found 
that adjusting the momentum grid removes the glitches, but we do not yet have a 
preventive fix. 
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Fig. 8. Deuteron wave functions in coordinate space for smooth and sharp cutoffs 
as in Fig. 7. 
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Fig. 9. The largest repulsive Weinberg eigenvalues at two energies as a function of 
cutoff for selected partial waves. The V\ ow k interactions are derived from the N 3 LO 
chiral potential of Ref. [16]. 

in the large-distance oscillations with increasing smoothing (for the S-state 
component this is visible only under additional magnification). Note that the 
"wound" in the N 3 LO coordinate-space wave function is removed by running 
down the cutoff (we plot u(r)/r to make the suppression near the origin more 
explicit). This feature of the wave functions leads to more perturbative be- 
havior in nuclear matter [5] as well as in few-body systems [8,10]. 

The "perturbativeness" of V\ ow k interactions with sharp cutoffs was exam- 
ined in Ref. [8] using Weinberg eigenvalues as a diagnostic. For details on the 
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Fig. 10. D~state probability Pd (left axis), binding energy Ed (lower right axis), 
and asymptotic D/S ratio rjd (upper right axis) of the deuteron as a function of the 
cutoff, starting from the Argonne vis [13] (left) and the N 3 LO chiral potential of 
Ref. [16] (right) with different smooth regulators. 

Weinberg analysis we refer the reader to Refs. [8,30]. In Fig. 9, we show the 
largest repulsive Weinberg eigenvalues as a function of the cutoff for selected 
channels, using the N 3 LO chiral potential from Ref. [16], which is constructed 
with a cutoff of 500 MeV. Although this is already a fairly soft potential, we 
still observe the characteristic decrease with low-momentum cutoff starting as 
high as 3.5 fm^ 1 (rather than at 2. 5fm~ 1 , as one might naively expect). This 
translates into weaker correlations in many-body wave functions (and there- 
fore better convergence). The rate of decrease is largely independent of the 
smoothness of the cutoff. 

Various deuteron properties are shown in Fig. 10 as a function of the cutoff 
using potentials derived from the Argonne fig potential [13] with exponential 
regulator n exp = 8 on the left, and from the N 3 LO chiral potential of Ref. [16] 
with exponential regulator n exp = 6 on the right. Plotted on the left axis is 
the D-state probability Pd, defined as 

roc rco 

P D = / drw{rf= / dkk 2 w(kf. (37) 
Jo Jo 

The cutoff dependence reflects the fact that Pd is not an observable [31,32]. 
The D-state probability evolves with the short-range part of the potential 
and, in particular, the short-range tensor interaction decreases as the cutoff is 
lowered. This decrease is desirable to reduce correlations in many-body wave 
functions [5]. The qualitative change in Pd is the same for other regulators 
and hermitization schemes (the Gram-Schmidt procedure is used in Fig. 10). 

In contrast to the D-state probability, the deuteron binding energy Ed and the 
asymptotic D/S ratio rjd are observables and thus cutoff independent (up to nu- 
merical tolerances), as shown by the right axes in Fig. 10. The asymptotic D/S 
ratio is calculated here by an extrapolation of the ratio of the deuteron wave 
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Fig. 11. The quadrupole moment of the deuteron calculated with the bare op- 
erator as a function of the cutoff for different smoothness regulators (left) and for 
different hermitization schemes (right). The low- momentum interactions are derived 
from the Argonne vis potential [13], and the experimental quadrupole moment is 
indicated with an arrow. 

function D to S-state components to the deuteron pole at k 2 = —mEd [31], 
rather than the conventional approach of extrapolating the 3 Si- 3 D! mixing 
angle ex. That is, we evaluate —w(k 2 )/u(k 2 ) on the Gauss mesh for positive 
k 2 , and then make a (near-linear) extrapolation to k 2 = —mE^. The constancy 
of rjd in Fig. 10 directly refutes the claim of Ref. [33] that this quantity cannot 
be reproduced using low-momentum interactions. 

Matrix elements of operators that are dominated by distance scales larger 
than the inverse cutoff are to a good approximation preserved as the cutoff is 
lowered. We investigate the evolution of operators by studying the expectation 
values in the deuteron of the bare quadrupole moment, rms radius and 1/r 
operators. The relevant formulas are [29] 



Qd = tzk drr 2 w{r) (V^w(r) - w{r)) 
20 Jo 
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Fig. 12. The rms radius of the deuteron calculated with the bare operator 
as a function of the cutoff for different regulators and hermitization schemes. 
The low-momentum interactions are derived using the N 3 LO chiral potential from 
Ref. [16]. 



and 



r°° 1 

(1/r) = J dr-[u(r) 2 + w(r) 2 ) 



(40) 



(We note that the momentum-space expressions for and show that these 
are only well-defined for a smooth cutoff [29].) These operators are dominated 
by the long-range part of the interaction, and therefore we expect that the 
expectation values change only for low cutoffs. This is verified in Fig. 11, where 
we show the deuteron quadrupole moment Qd as a function of the cutoff for 
different smoothness regulators (left) and for different hermitization schemes 
(right), and in Figs. 12 and 13, where the rms radius and the matrix element 
of 1 jr are plotted ClS db function of A. 



We emphasize that the bare quadrupole moment by itself does not correspond 
to an experimental observable. To correctly reproduce the experimental mo- 
ment (indicated in Fig. 11 by an arrow), one needs the corresponding op- 
erators. We find that the change of the quadrupole moment is of the same 
magnitude as the difference between the bare and experimental quadrupole 
moments for cutoffs down to 1.5-2 fm _1 . From the deuteron wave functions 
shown in Fig. 14, we conclude that low-energy observables in the deuteron 
are reproduced without short-range correlations in the wave function. Similar 
observations hold for the expectation values of the rms radius and of the 
1/r operator. 

The key motivation for smooth cutoffs was to remedy the slow convergence 
at the lOkeV level in the deuteron and at the 100 keV level in the triton, 
when calculated in a harmonic oscillator basis. In Fig. 15, the relative error 
in the binding energy of the deuteron (with respect to the converged result) 
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Fig. 13. The matrix element of 1/r in the deuteron calculated with the bare oper- 
ator as a function of the cutoff for different regulators and hermitization schemes. 
The low-momentum interactions are derived using the N 3 LO chiral potential from 
Ref. [16]. 
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Fig. 14. Deuteron wave functions in coordinate space for low-momentum interactions 
at several different cutoffs using an exponential regulator with n exp = 8. The initial 
interaction is the N 3 LO chiral potential from Ref. [16]. 

is shown as a function of the size of the oscillator space (iV max hw excitations) 
for a range of regulators. The slow convergence is evident for the sharp cutoff, 
where the error is below the percent level only for the largest space. The Fermi- 
Dirac regulator with e = 0.2 fm" 1 , which is still very sharp (see Fig. 2), gives 
improved errors but still requires large spaces. Increasing e from 0.2 fm -1 to 
0.5 fm -1 steadily decreases the error until it improves rapidly with the space 
size, while still only minimally distorting phase shifts. Very similar relative 
errors are obtained for the exponential regulator with n = 8 (not shown) or 
n = 10. Lowering n to 4 reduces the error still further, but at the cost of a 
potentially significant distortion of the phase shifts. 
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Fig. 15. The relative error in the deuteron binding energy as a function of the 
size of the oscillator space for sharp cutoff and various smooth regulators. 
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Fig. 16. The absolute error in the deuteron binding energy Ed calculated from 
low-momentum interactions derived using the Argonne v\% potential [13] (left) and 
the N 3 LO chiral potential from Ref. [16] (right) as a function of the cutoff. We show 
results for two oscillator spaces and for two different exponential regulators. 

The absolute error in the deuteron binding energy is shown in Fig. 16 as a 
function of the cutoff for low-momentum interactions derived from the Ar- 
gonne t>ig [13] and the N 3 LO chiral potential from Ref. [16]. In all cases, the 
reproduction of the binding energy improves with decreasing cutoff until the 
limits of numerical accuracy are reached. This is in contrast to the observed 
degradation below 2 fm -1 for a sharp cutoff [9]. The improvement for the N 3 LO 
potential starts just below 4fm _1 (even though the EFT cutoff is 500 MeV or 
2.5 fm~ ) and is quite significant by 2fm~ . 



23 




The convergence is also greatly improved for the triton. In Fig. 17, the triton 
binding energy is plotted as a function of the size of the oscillator space. For 
efficiency, convergence for the smallest possible space is desirable, as the com- 
putational cost grows rapidly with iV max and for larger systems. Convergence 
at the keV level is achieved by all exponential regulators with n ^ 4 (n exp = 8 
is shown) soon after iV max = 20. The consequence in moving from sharp to 
increasingly smooth regulators is seen from the Fermi-Dirac regulators, where 
e = 0.5 fm" 1 yields very satisfactory results. The difference between converged 
results for the triton is a measure of the differences in the short-range three- 
body force with the different regulators. 

The dramatically improved convergence of lower cutoffs compared to even the 
soft N 3 LO chiral potential is demonstrated in Fig. 18. Note that the phase 
shifts would be essentially undistorted at A = 2.5 fm . In addition, when we 
evolve the potential to low cutoffs, the dependence on the oscillator parameter 
b becomes flatter for a given iV max , and for A = 2.0 fm -1 the minimum is 
unchanged as the space is enlarged. The sharp-cutoff V\ ow k also shows a flatter 
dependence on the oscillator parameter compared to the initial N 3 LO chiral 
potential and converges rapidly to the 100 keV level, but then converges slowly 
(the convergence for the N 3 LO chiral potential is very similar to that shown 
for the Argonne Vis potential in Fig. 17). 

In Fig. 19, the triton binding energy as a function of the cutoff is shown with 
the exponential regulator n exp = 6 calculated from the Argonne Vi$ [13] and 
using the N 3 LO chiral potential of Ref. [16]. For comparison, we also show the 
Faddeev results from Ref. [4], which use a sharp cutoff. The triton binding 
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oscillator parameter b [fm] 

Fig. 18. Dependence of the triton binding energy E% as a function of oscillator 
parameter b for the bare N 3 LO chiral potential of Ref. [16] (upper left) and evolved 
to low momenta A = 3.0,2.5 and 2.0 fm -1 with exponential regulator n exp = 8. 
Results are presented for three smaller spaces -/V max = 8, 12 and 16, along with the 
converged iV max = 40 energy. 

energy is cutoff dependent because we have not included the corresponding 
three-body interactions. Therefore, the difference to the experimental binding 
energy (shown with an arrow) is the total three-body contribution. The pat- 
tern of running is very similar in the three cases, although the values for a 
given regulator naturally vary. While it is possible to choose a cutoff so that 
the experimental triton binding energy is reproduced by only the two-body 
interaction, the total three-body contribution will not vanish in other systems 
(e.g., the alpha particle [4] or nuclear matter [5]). 



5 Summary and Conclusions 

In this paper, low- momentum "V[ ow fc" interactions with smooth cutoffs were 
constructed and tested for several different types of regulators and hermiti- 
zation schemes. Problems seen in previous work with variational calculations 
that were attributed to artifacts from a sharp cutoff are resolved by smooth 
regulators. In particular, convergence for the energy of the deuteron and tri- 
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Fig. 19. The triton binding energy E t as a function of the cutoff starting from the 
Argonne v\% potential [13] with a sharp cutoff (taken from the Faddeev calculations 
of Ref. [4]) and with an exponential regulator. We also show results for the N 3 LO 
chiral potential of Ref. [16] with the same exponential regulator. 

ton in a harmonic oscillator basis is greatly improved. We have checked that 
this conclusion also holds for other variational trial wave functions not shown 
here. We expect this improved convergence to carry over to other few- and 
many-body systems. 

The regulators introduced here are specified not only by the cutoff on relative 
momenta but by a parameter that determines the smoothness. Our conclu- 
sions for the optimal regulator are provisional, because this work has been 
restricted to convergence rates in the deuteron and the triton. However, based 
on the observed convergence it appears that a similar degree of smoothness is 
optimal for different regulators. This is achieved with n ~ 8 for the exponen- 
tial regulator or e ~ 0.5 fm _1 for the Fermi-Dirac regulator. Going to smoother 
regulators maintains the same rapid convergence but will increase distortions 
of the phase shifts. Whether or not this is an issue depends on the application, 
but we expect that the latter is not a problem for low-energy observables. 

Given the regulator freedom that exists in generating low-momentum interac- 
tions for a given cutoff A, it is clearly a misnomer to speak of "the" V\ ow k poten- 
tial. We emphasize in particular that the consistent three-body (and higher 
many-body) interactions corresponding to the smooth two-body interaction 
will differ in each case. More precisely, since the long-distance (pion-exchange) 
parts of the interaction are preserved by the RG (until the cutoff is compa- 
rable to the pion mass), the short-distance part of the three-body interaction 
is modified. Cutoff and regulator independence with two-body interactions 
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alone should not be a criterion for the applicability of low-momentum inter- 
actions to nuclear structure. (Note that this perspective differs markedly from 
that expressed in Ref . [34] , which discussed the application of low- momentum 
potentials with sharp cutoffs to few-nucleon systems.) Rather, the three-body 
variations for different regulators allow us to extend the powerful test of cutoff 
independence for observables to include independence of the regulator. How- 
ever, without a consistent three-body interaction, few- and many-body cal- 
culations are largely meaningless (except for considerations of convergence). 
Therefore, a high priority for the near term is to develop the machinery to 
quickly fit approximate three-body interactions for a given regulator. 

The methods described here apply equally to low-momentum interactions de- 
rived from conventional nucleon-nucleon potentials and to those derived from 
chiral EFT potentials. The latter has the advantage of a systematic organiza- 
tion of many-body forces and operators. While chiral EFT potentials already 
start from lower cutoffs than conventional NN interactions, we find significant 
added advantages for few- and many-body calculations by starting with chiral 
potentials fit at a larger cutoff and running them down to a lower cutoff [5,8]. 

Field redefinitions in EFT reshuffle higher-order terms so that the truncation 
error is different but of the same order. That is, if the EFT is specified to NLO, 
then N 2 LO terms will differ after a field redefinition but the expected trunca- 
tion error is still N 2 LO. The RG transformations discussed here preserve the 
error as the cutoff is lowered by generating all necessary higher-order short- 
range operators through the evolution. This is particularly advantageous when 
the running is rapid and the separation of scales is not large, as in the nuclear 
case (when the tensor force is active). At present, only the two-body interac- 
tion is evolved while the three-body interaction is fit for each cutoff using the 
N 2 LO form [4]. While there are solid indications that this is a reasonable pro- 
cedure [4], the evolution of consistent chiral three-body interactions to lower 
momenta is a major goal. In addition, the smooth regulators discussed here 
may have advantages for constructing chiral EFT potentials with low cut- 
offs. As is apparent from Fig. 2, the exponential regulator with n 8 and 
the Fermi-Dirac regulator with e ~ 0.5 fm -1 lead to weaker distortions than 
the conventional exponential regulator with n = 3 in the N 3 LO chiral poten- 
tials [16,17], and thus smaller errors have to be absorbed by the counterterms. 

There are immediate applications that can be made using the low-momentum 
interactions with smooth cutoffs developed here. Variational calculations using 
hyperspherical harmonics [35], the No-Core Shell Model [36] and the Coupled- 
Cluster approach [37] should benefit from the improved convergence. Tests for 
all of these are in progress. 
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